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•'— I ' Abstract 

1 " I ' We derive equations of motion for the dynamics of anisotropic particles directly from the dis- 

sipative Vlasov kinetic equations, with the dissipation given by the double bracket approach 
' (Double Bracket Vlasov, or DBV). The moments of the DBV equation lead to a nonlocal form 

of Darcy's law for the mass density. Next, kinetic equations for particles with anisotropic 
interaction arc considered and also cast into the DBV form. The moment dynamics for 
, these double bracket kinetic equations is expressed as Lie-Darcy continuum equations for 

' densities of mass and orientation. Wc also show how to obtain a Smoluchowski model from a 

— I . cold plasma-like moment closure of DBV. Thus, the double bracket kinetic framework serves 

' as a unifying method for deriving different types of dynamics, from density-orientation to 

. Smoluchowski equations. Extensions for more general physical systems are also discussed. 

o ■ 

> ■ 

■ 1 Introduction 

X: 

■ 1.1 Geometric models of dissipation in physical systems 

This paper explains how the geometry of double-bracket dissipation makes its way from the 
microscopic (kinetic theory) level to the macroscopic (continuum) level in the process of taking 
moments of the Vlasov probability distribution function, when the particles in the microscopic 
description carry an internal variable that is orientation dependent. Without orientation de- 
pendence, the moment equations derived here yield a nonlocal variant of the famous Darcy law 
|Darcyl856| . When orientation is included, the resulting Lie-Darcy moment equations identify 
the macroscopic parameters of the continuum description and govern their evolution. 

In previous work. Gibbons, Holm and Kupershmidt |GiHoKul982l IGiHoKul983j (abbrevi- 
ated GHK) showed that the process of taking moments of Vlasov's equation for such particles 
is a Poisson map. GHK used this property to derive the equations of chromohydrodynamics. 
These are the equations of a fluid plasma consisting of particles carrying Yang-Mills charges 
and interacting self-consistently via a Yang-Mills field. The GHK Poisson map for chromohy- 
drodynamics was an extension of the Kupershmidt-Manin (KM) bracket |KuMal978j for the 
moments in the integrable hierarchy of long- wave equations due to Benney |Bel973j . Gibbons 
|Gil981j had noticed that the KM Poisson map for the Benney equations would also apply to 
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the moments of a Vlasov (collisionless) plasma |Vlasovl961j for particles carrying an ordinary 
electric charge. GHK extended the KM bracket to allow the particles to carry orientation (or 
Yang-Mills) properties. GHK considered only Hamiltonian motion. They did not consider the 
corresponding double-bracket Poisson structure of dissipation. That is the subject of the present 
work. 

1.1.1 History of double-bracket dissipation 

Bloch, Krishnaprasad, Marsden and Ratiu ( |BlKrMaRal996] abbreviated BKMR) observed that 
linear dissipative terms of the standard Rayleigh dissipation type are inappropriate for dynam- 
ical systems undergoing coadjoint motion. Such systems are expressed on the duals of Lie 
algebras and they commonly arise from variational principles defined on tangent spaces of Lie 
groups. A well known example of coadjoint motion is provided by Euler's equations for an ideal 
incompressible fluid |Arl966j . Not unexpectedly, adding linear viscous dissipation to create the 
Navier-Stokes equations breaks the coadjoint nature of the ideal flow. Of course, ordinary vis- 
cosity does not suffice to describe dissipation in the presence of orientation-dependent particle 
interactions. 

Restriction to coadjoint orbits requires nonlinear dissipation, whose gradient structure differs 
from the Rayleigh dissipation approach leading to Navier-Stokes viscosity. As a familiar example 
on which to build their paradigm, BKMR emphasized a form of energy dissipation (Gilbert 
dissipation |Gilbertl955] ) arising in models of ferromagnetic spin systems that preserves the 
magnitude of angular momentum. In the context of Euler-Poincare or Lie- Poisson systems, this 
means that coadjoint orbits remain invariant, but the energy decreases along the orbits. BKMR 
discovered that their geometric construction of the nonlinear dissipative terms summoned the 
double bracket equation of Brockett [Brl988t IBrl993j . In fact, the double bracket form is well 
adapted to the study of dissipative motion on Lie groups since it was originally constructed as 
a gradient system |Blochl99n[ IBlBrR,al992l IBr 19941 [BlFlR,al996j . 

While a single Poisson bracket operation is bilinear and antisymmetric, a double bracket 
operation is a symmetric operation. Symmetric brackets for dissipative systems, particularly for 
fluids and plasmas, were considered previously by Kaufman |Kal9841 IKal985) . Grmela |Grl9841 
I(kl993a[ IGrl993bj . Morrison |Mol984[ IMol986j . and Turski and Kaufman |TuKal987j . The 
dissipative brackets introduced in BKMR were particularly motivated by the double bracket 
operations introduced in Vallis, Carnevale, and Young [VaCaYol989] for incompressible fluid 
flows. 

1.1.2 Selective decay hypothesis 

One of the motivations for Vallis et al. |VaCaYol989j was the selective decay hypothesis, 
which arose in turbulence research |MaMol980| and is consistent with the preservation of coad- 
joint orbits. According to the selective decay hypothesis, energy in strongly nonequilibrium 
statistical systems tends to decay much faster than certain other ideally conserved properties. 
In particular, energy decays much faster in such systems than those "kinematic" or "geometric" 
properties that would have been preserved in the ideal nondissipative limit independently of the 
choice of the Hamiltonian. Examples are the Casimir functions for the Lie-Poisson formulations 
of various ideal fluid models |IIoMaRaWel985] . 

The selective decay hypothesis was inspired by a famous example; namely, that enstrophy 
decays much more slowly than kinetic energy in 2D incompressible fluid turbulence. Kraichnan 
|Krl967j showed that the decay of kinetic energy under the preservation of enstrophy causes 
dramatic effects in 2D turbulence. Namely, it causes the well known "inverse cascade" of kinetic 
energy to larger scales, rather than the usual "forward cascade" of energy to smaller scales 
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that is observed in 3D turbulence. In 2D ideal incompressible fluid flow the enstrophy (the 
norm of the vorticity) is preserved on coadjoint orbits. That is, enstrophy is a Casimir of 
the Lie-Poisson bracket in the Hamiltonian formulation of the 2D Euler fluid equations. Vallis 
et al. |VaCaYol989j chose a form of dissipation that was expressible as a double Lie-Poisson 
bracket. This choice of dissipation preserved the enstrophy and thereby enforced the selective 
decay hypothesis for all 2D incompressible fluid solutions, laminar as well as turbulent. 

Once its dramatic effects were recognized in 2D turbulence, selective decay was posited as a 
governing mechanism in other systems, particularly in statistical behavior of fluid systems with 
high variability. For example, the slow decay of magnetic helicity was popularly invoked as a 
possible means of obtaining magnetically confined plasmas |Tal986j . Likewise, in geophysical 
fluid flows, the slow decay of potential vorticity (PV) relative to kinetic energy strongly influences 
the dynamics of weather and climate patterns much as in the inverse cascade tendency in 2D 
turbulence. The use of selective decay ideas for PV thinking in meteorology and atmospheric 
science has become standard practice since the fundamental work in |HoMcRol985t IYol987| . 

1.1.3 Kandrup and the double bracket for astrophysical instabilities 

A form of selective decay based on double-bracket dissipation was proposed in astrophysics by 
Kandrup |Kal99lj for the purpose of modeling gravitational radiation of energy in stars. In 
this case, the double-bracket dissipation produced rapidly growing instabilities that again had 
dramatic effects on the solution. The form of double-bracket dissipation proposed in Kandrup 
[Kal99lj is a strong motivation for the present work and it also played a central role in the 
study of instabilities in BKMR. 

The double bracket idea mentioned earlier in the context of magnetization dynamics |Gilbertl955j 
implements dissipation by the sequential application of two Poisson bracket operations. The dis- 
sipative Vlasov equation in |Kal99lj is written as: 



where f{q,p, t) is the probability distribution function for a single particle in phase-space, a > 
is a positive constant, H is the Vlasov Hamiltonian (depending on the system under consider- 
ation) and {• , •} is the canonical Poisson bracket in {q,p)- Here we use standard notation in 
kinetic theory: in the general case, q is a point on a Riemannian manifold Q, while {q,p) ^ T*Q 
is a point in the cotangent bundle. In common situations one considers Q = M'^, so that 
T*Q = M^'^. In this paper we shall consider a one dimensional physical space (Q = M) to 
simplify the treatment. 

When a — 7> in (jl.ip . this equation reduces the Vlasov equation for collisionless plasmas. 
For a > 0, this is the double bracket dissipation approach for the Vlasov-Poisson equation. 
This nonlinear double bracket approach for introducing dissipation into the Vlasov equation 
differs from the standard Fokker-Planck linear diffusive approach |Fokker-Plankl93r| . which 
adds dissipation on the right hand side of equation (jl.ip as a linear term given by the Laplace 
operator in the momentum coordinate Apf. 

1.1.4 The double bracket and Riemannian geometry 

An interesting feature of the double bracket formulation is that it leads via a variational approach 
to a symmetric Leibnitz bracket that in turn yields a metric tensor and an associated Riemannian 
(rather than symplectic) geometry for the solutions. The variational approach thus preserves 
the nature of the evolution of Vlasov phase space density, by coadjoint motion under the action 
of the canonical transformations on phase space densities. 
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As Otto |Ot2001j explained, the Riemannian geometry of dissipation may be revealed by un- 
derstanding how it emerges from a variational principle. Here, we follow the variational approach 
and consider a generalization of the double bracket structure in equation (jl.ip that recovers pre- 
vious cases for particular choices of modeling quantities introduced in [HoPuTr2007-CR] . 



The important feature of this double-bracket approach to dissipation (which is necessarily 
nonlinear) for the Vlasov equation is that this form of the dissipation preserves the most ba- 
sic property of Vlasov dynamics. Namely, the dissipation is introduced as a canonical 
transformation. This means that the solution is still an invariant probability distribu- 
tion and satisfies Liouville's theorem, even though its dynamics do not conserve either 
the Hamiltonian H or the energy E. Eq. (II. 2p generalizes the double bracket operation in 
Eq. (jl.ip and reduces to it when the Hamiltonian H is identical to the dissipated energy E 
and the mobility fi[f] = af is proportional to the Vlasov distribution function / for a positive 
constant q > 0. As shown in jHoPuTr2007-CR] the generalization ()1.2p has important effects 
on the types of solutions that are available to this equation. Indeed, the form (II. 2p of the 
Vlasov equation with dissipation allows for more general mobilities than those considered in 
|BlKrMaRal9961 IKal99l] IKal984l IMol984j . For example, one may choose mobility in the form 
fj,[f] = K*f, where the ( * ) operation is a convolution in phase space with an appropriate kernel 
K. In [HoPuTr2007j the smoothing operation in the definition of p[f] introduces a fundamental 
action scale (the area in phase space associated with the kernel K) into the dissipation mech- 
anism. Accordingly, the dissipation depends on phase-space averaged quantities, rather than 
local pointwise values. 

This smoothing also has the fundamental advantage of endowing (jl.2p with the one-particle 
solution as its singular solution. The generalization Eq. (jl.2p may also be justified by using 
thermodynamic and geometric arguments |HoPuTr2007j . Section [2] shows that this generaliza- 
tion leads from the microscopic kinetic level to the classic Darcy law (velocity being proportional 
to force) in the continuum description. 

1.2 Our goal and approach 

The goal of the present work is to determine the macroscopic implications of introducing non- 
linear double-bracket dissipation at the microscopic level, so as to respect the coadjoint orbits 
of canonical transformations for dynamics that depends upon particle orientation. Our ap- 
proach introduces this orientation dependence into the microscopic description by augmenting 
the canonical Poisson bracket in position q and momentum p so as to include the Lie-Poisson 
part for orientation g taking values in the dual g* of the Lie algebra q, with eventually q = so (3) 
for physical orientation. Here we denote by Q the Lie group underlying the particle symmetry 
properties, while the corresponding Lie algebra g = T^Q and its dual g* = T*Q are denoted by 
gothic fonts. Although it is unusual to denote a point in g* by g, we shall keep this notation to 
avoid confusion with the Greek indexes used below. 

The identification of the particle orientation g with an element in 50* (3) is consistent with 
ordinary theories in condensed matter physics, such as the theory of sin glasses. In this context, 
the spin angular momentum variable s £ 50* (3) possesses the typical Poisson bracket {s^, Si,} = 
^iiv\S\ |DzVol980] . where the Levi-Civita symbol e denotes the total antisymmetric tensor. 
In more general situations when the rotational symmetry is completely broken, one also needs 
to consider the particle rotational state G € 5*0(3) as an extra Lie group coordinate in the 
augmented phase-space. This is relevant, for example, when one wants to consider frustration 
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effects. However, here we shah not deal with symmetry breaking and we shall simply assume that 
the orientational degree of freedom can be identified with the momentum variable g € 50* (3), 
or in more generality g € Q* . 

Prom the considerations above, one concludes that the introduction of the orientational 
degree of freedom g leads to a Vlasov distribution f{q,p,g,t) on the augmented phase space 
T*Q X 5*, whose total Poisson bracket is written in the GHK form as 



f,h]^ := [f,h] 



+ {g, 



(1.3) 



"9/ dh' 

dg' dg 

where [■, •] : gxg— )-0is the Lie algebra bracket and ( • , • ) : g* x g — )■ R is the pairing between 
the Lie algebra g and its dual g*. For rotations, g = so(3) and the bracket [• , •] becomes the 
cross product of vectors in R^. Correspondingly, the pairing ( • , • ) becomes the dot product of 
vectors in R^. We consider the double-bracket dynamics of f{q,p,g,t) resulting from replacing 
the canonical Poisson brackets in Eq. (II. 2p by the direct sum of canonical and Lie-Poisson 
brackets { • , • }i in Eq. (|1.3p . We then take moments of the resulting dynamics of f{q,p,g,t) 
with respect to momentum p and orientation g, to obtain the dynamics of the macroscopic 
description. The moments with respect to momentum p alone provide an intermediate set of 
dynamical equations for the p— moments, 

An{q,g,t) := J p"" f{q,p,g,t)dp with ^ G g* . (L4) 

(In higher dimensions one replaces the power with the tensor power p®", also denoted by 
(8>"p.) These intermediate dynamics are reminiscent of the Smoluchowski equation for the prob- 
ability ^o('i') 9) t)- We shall show how the derivation of usual Smoluchowski-type equations from 
double-bracket arguments may require special care depending on the chosen moment closure. 
In particular, we formulate a new Smoluchowski equation, which encodes the Landau-Lifshitz- 
Gilbert dissipation for ferromagnetic particles. 

We also find a closed set of continuum equations for the (p, g) moments of the type 

AnAQ^t) := / p^'g-' f{q,p,g,t)dpdg. (L5) 



where g'^ denotes the tensor power g'^ = g ® ■ ■ ■ ® g, u times. The final closure provides the 
macroscopic continuum dynamics for the set of moments of the double-bracket kinetic equations 
(jl.2p under the replacement {•, •}—>■{•, - ji. This macroscopic continuum closure inherits the 
geometric properties of the double bracket, because the process of taking these moments is a 
Poisson map, as observed in GHK |GiHoKul982l IGiHoKul983l IGiHoTr2008] . 



1.3 Mathematical framework for geometric dissipation 

As explained in |Ot2001j dissipation of energy E may naturally summon an appropriate metric 
tensor. In previous work Holm and Putkaradze |HoPu20071 IHoPuTr2007j introduced a distance 
between functions on a vector space V that is acted on by a Lie algebra X of vector fields. 
In particular they showed that that for any two functionals of a geometric tensor 

quantity k G y a distance between them may be defined via the Riemannian metric. 




Here ( • , • ) denotes the pairing of vector fields in X with their dual one-form densities in 
X*, sharp (•)" raises the vector index from covariant to contravariant and //[k] is the mobility 
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functional. The mobility h[k\ is assumed to satisfy the requirements for (|1.6|) to be positive 
definite and symmetric, as discussed in |HoPuTr2007j . The diamond operation (o) in equation 
()1.6p is the dual of the Lie algebra action, defined as follows. Let a vector field ^ € X act on a 
vector space V by Lie derivation, so that the Lie algebra action of ^ on any element k (z V is 
given by the Lie derivative, ^- k = £^k. The operation dual to the Lie derivative is denoted by 
o and defined in terms of the pairings between X and X* and between V and V* as 

(CoK,d :=(c,-4^) ■ (1-7) 

Given the metric (|1.6p and a dissipated energy functional E[k], the time evolution of arbitrary 
functional F[k] is given by |HoPu2007[ IHoPuTr2007] as 

F = {{F,E}}[K]:=-g,{F,E) , (1.8) 

which specifies the dynamics of any functional F[k], given the the energy dependence E[k]. The 
bracket {{F , E}} is shown to satisfy the Leibnitz product-rule property for a suitable class 
of mobility functional //[k] in |HoPu2007l lHoPul'r2007] . Eq. ([TSD and positivity of gn{E,E) 
imply that the energy E decays in time until it eventually reaches a critical point, 5E/5k = 0. 
For densities (dual to functions in the pairing), the Lie derivative is the divergence: £^k = 
div(K,^) and its dual operation is (minus) the gradient: k o = — kVC- Thus, for densities the 
symbol diamond ( o ) is replaced by gradient ( V ) in the metric defined in (II. Sh . 

The definition of the dissipative bracket in Eq. (jl.Sp for arbitrary functionals {{F , E}} [k] is 
the basis for our present considerations of dissipation in kinetic equations. In this paper we will 
extend the geometric dissipation (jl.Sp to the symplectic case by defining the star {*) operator, 
which is the analogue of diamond (o) for symplectic spaces. 



2 Dissipation in the Vlasov equation and the moment hierarchy 
2.1 Dissipative bracket for the Vlasov equation 

The dissipative term in equation (II. 2p is found by considering the symplectic algebra action of 
a Hamiltonian vector field associated with a Hamiltonian function h. This action is given 
on a phase space density / through the canonical Poisson bracket {■ , •} as follows |MaWe82] : 
f = {/) ^} =■ h ■ f. One can check that the dual diamond operator in ()1.7p (denoted by 
■k for the symplectic Lie algebra action above) is still a Poisson bracket |HoPuTr2007j . so that 
f = {g, /}. Thus, we introduce the dissipative bracket (11.81) 

«.,.H^-{{./l,f}',{/,e)^({/,{./l,f}'},^) ... 

which gives the dissipative term in equation (jl.2p . Here, the sharp ft operator is a suitably chosen 
mapping from the (dual) space of distributions to the Lie algebra of phase space functions, so 
the pairing can be taken. Usually, we take the pairing to be a simple pairing, in which case 
the sharp notation is unnecessary, as in the previously postulated form of the equations (|1.2p . 
The sharp notation in (12. 1|) is used to indicate that the pairing may involve a metric. 

Assumption 2.1 Everywhere in this paper, we assume that the Lie algebra q of orientation 
motions is finite dimensional. Of particular interest to application is the motion of magnetized 
particles, described by SO{3) group. For these motions, q = so(3) ~ M'^, so the integrals with 
respect to g (z Q are simply three-dimensional integrals. 
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Extensions of our theory could be computed for an arbitrary finite- dimensional Lie algebras, 
describing, for example, particles with a finite number of joints. However, in this paper we will 
stay away from infinite-dimensional Lie algebras, that may arise, for example, for the case of 
flexible filaments. Even though the formal extension of our double bracket dissipation theory is 
possible even in that case, the mathematical difficulties of defining integrals and corresponding 
measures are beyond the scope of this paper. On the other hand, the canonical variables mo- 
menta p and coordinates q describe the motion of the center of mass of the particles and thus 
the integrals with respect to p and q are always finite-dimensional. In order to elaborate further, 
we introduce the following 

Assumption 2.2 Everywhere in the paper we assume that all the finite-dimensional integrals 
defining the moments U.5\) exist. Formally, this can be achieved, for example, by requiring that 
the distribution function f is smooth and decays to zero sufficiently fast whenever p, q or g tend 
to oo in some properly defined norm. While for compact configuration manifolds requirements 
on q integration convergence are automatically satisfied, the convergence requirements still need 
to be enforced for p and g. 

An interesting point about this form of double-bracket dissipation in kinetic equations is 
that it leads rather naturally to another widely used dissipative equation; namely, the venerable 
Darcy law |Darcyl856| . In particular, we will show in this section how one recovers the nonlocal 
version of Darcy's law introduced in |HoPu2005t IHoPu2006] . This could be accomplished by 
integrating the purely dissipative kinetic equation 

f ^ {/,{.!/,, f}}, 

with respect to the momentum coordinate p (here we have dropped the sharp notation for 
simplicity). We shall call this equation Double Bracket Vlasov equation (DBV). Note however, 
that this DBV equation is valid for isotropic particles only. More general version of DBV 
equation including the anisotropy will be derived below in (j3.ip . In order to derive Darcy's law 
for particle motion, we apply the geometric dissipative bracket (|1.8|) to the equations of the 
moment hierarchy. 

2.2 The Kupershmidt-Manin bracket for ID moment dynamics 

As a general result |Gil981j . the equations for the moments of the Vlasov dynamics (Eq. 
(jl.ip with a —7- 0) form a Lie-Poisson system under the Kupershmidt-Manin (KM) bracket 
[KuMal978j . The n-th moment is defined as An{q, t) := J f{q,p, t) dp. These quantities have 
a geometric interpretation in terms of (covariant) tensor densities |GiHoTr2005t IGiHoTr20071 
IGiHoTr2008j which can be seen by re-writing the moments as An = &^{pdq) f{q,p) dq A dp= 
An{q) ^^"^^dq, where 'S'^'dq := dq ^ ■ ■ ■ dq n times. Thus, moments An belong to the vector 
space dual to the contravariant tensors of the type /3n = fin{(l) <8)"i9g. Upon denoting the spatial 
derivative d/dq by prime, these tensors are given the following Lie algebra structure |Nil955| 

{am, /?„! = ( n/3„(<?) «:„((?) -ma^(g)/3^(<7)) =: ad„„ /3„ (2.3) 

so that the Kupershmidt-Manin Poisson bracket for moment dynamics is |KuMal978tlGiHoTr2008"] 

{F^H] = l^Am+n-l, 

where we sum over repeated indices. Thus, for a given moment Hamiltonian H, the Vlasov 
moment equations are dtAn = — ad|j:^/^^^ ^m.+n-i, where the ad* operator is defined by 



5F 5H 

6An SAni 



(2.4) 
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( ad^^ Ak, ak-n+i ) ■= {Ak, adp„ ak-n+i ) and is given explicitly as ad^^ = ((A; + 1) A^ dg /3„ + 
n/3ndqAk) (8''^""''''^ dq. In the next section we use the following relation between the KM ad* 
operator and the canonical Poisson bracket 

m=0 ^ ^ ■> m=0 

where we have used the chain rule formula 5H /5f = p™ 6H /6Am, that is ordinarily used to pass 
from the Vlasov description to its corresponding moment formulation |Gil981t IGiHoTr2008] . 
For completeness, we report the following basic fact underlying the Hamiltonian structure of ID 
moment dynamics 

Theorem 2.3 (Gibbons |Gil981| ) The operation of taking moments of Vlasov solutions is a 
Poisson map. 

This theorem was recently extended to any number of dimensions and to anisotropic interactions 
in [GiHoTr2008j . where the Poisson map property is justified by the use of momentum maps. 

2.3 Geometric dissipation for moment dynamics and Darcy's Law 

Consider the following Lie algebra action on Vlasov densities: /3„ • / := £xpn^^/ = {/> P^Pn} 
(no sum), which is naturally given by the canonical action of the Hamiltonian function h{q,p) = 
p"'f3n{q)- The dual action defines the -kn operator, given by 

{f^ng, Pn) := (/, Pn- g) = j f {9 , p"" Pn} dq dp = (^j {f,g}p''dp, /3„^ . (2.5) 

Consequently, the dissipative bracket for the moments is written as |HoPuTr2007-CR] 

dF / f „ ( 6E^ , f „ ( „ SF' 

dt 



{{E,F}} = -Mp«|;.[/],^|dp, y"p"|/,^|dp 



Upon writing SE/Sf = p^Pk and SF/Sf = p^am, the dissipative bracket becomes 

{{E,F}} = -(^d*p^]lk+n-l, (ad:^^^+„_i)«) , (2.6) 
where fi-sio) '■= Jp^t^lf] dp. The purely dissipative dynamics for the moments is then given by 

= ad*^An+n-i with 7„ := (ad^^, /Ifc+m-i)" . (2.7) 

Here, the 7^ are the tensor fields representing the 'velocities' for the tensor densities Am+n-i- 
It turns out that the equation for Aq carried along by the velocity vector field 71 is exactly 
Darcy's law. Indeed, upon dropping the sharp for simplicity and by truncating equation (|2.7|) 
to consider only m = 1, we obtain the equation dtp = -^71 p = dq{p^i), for the density p := Aq. 
If we assume that E = E[p\ to neglect inertial effect and also that p^ = p[p\ for simplicity, then 
71 = (podPo/dq)'^ is recognised as the Darcy velocity. Therefore, the density p := Aq evolves 
according to Darcy's law, namely 

dp d f d 6E\ 

m=d-qV^^P^d-qJ^) ■ ^'-'^ 

Notice that no linear diffusion terms are present in the above equation, since such terms can 
only arise from Brownian particle motion that is not considered in this paper. Linear diffusion 
terms can however be added a posteriori to account for stochastic effects. 
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Remark 2.4 (Local vs. nonlocal mobilities) We notice that in all our discussions it is not 
necessary for the mobility quantity fj, to be a functional of the dynamical variable. Indeed, even in 
the general formulation (II. 6|) of double-bracket dissipation, the quantity //[k] must only satisfy 
the requirements for (jl.6p to be symmetric and positive definite. Thus, at the kinetic level, 
the phase-space mobility //[/] in the dissipative Vlasov equation ()2.2p can also be an ordinary 
function fi^f) (rather than a functional nlf]) of the Vlasov distribution f{q,p) or even some 
fixed phase-space density fi{q,p). At the macroscopic level of Darcy's law, this would correspond 
respectively either to a function fi{p) of the macroscopic particle density p{q) or to some fixed 
function p{q) on physical space. In these cases, the double bracket Vlasov equation (j2.2p simply 
produces Darcy's Law with complete absence of nonlocal terms. However, in what follows we 
shall consider the more general choice of a mobility functional, whose corresponding nonlocal 
effects have been investigated in JHoPu2005^ . 

2.4 Vlasov derivation of dissipative moment equations 

This section is devoted to clarify how one can pass from the dissipative Vlasov equation (j2.2p 
to the double-bracket equations for the moments ()2.7p . The starting point is to recognize that 
the double-bracket structure of equation (j2.2p involves a sharp operation {■)^ as follows 

0.1.1/1,^1 VI. (2..) 



dt y ' r 

This sharp operation is evidently defined through an appropriate (co-)metric, and thus the most 
general form of the double-bracket Vlasov equation ()2.2p is written as 



dt ' i"^^^' 6f 

Here K = K{q,p,q' ,p') is a symmetric positive-definite kernel, so that the sharp operator is 
given by the convolution K * f = f^ taking a density variable / into a scalar function The 
simplest choice of a 5-function for K evidently returns equation ()2.2p . In a more general case, 
one can consider a kernel K that is analytic in the momentum coordinates, so that the double 
bracket structure 



5f \ i Sf ] / JJi ^f\\ V 5f 



(2.11) 



becomes 



where we have used the Taylor expansion 

oo 

K{q,p,q',p')= J2 p''p'"'Kr.,m{q,q')- (2-13) 



n,m=0 



Assumption 2.5 In what follows, we shall assume that the coefficients Kn,m decay sufficiently 
fast so that the convergence of all the necessary sums is guaranteed. 
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Upon expanding 5F/6f = "^p^ 5F / 5Ai (and analogously for 6E/6f) and integrating over the 
momentum variables, one obtains the double-bracket structure for the moments in the form 

{{E,F}} = - ^ (ad;^ Jlk+n-i , Kn,i * {adl^Am+i^i) > , 

n,l,k,m 



which evidently recovers (|2.6p . 

Upon setting k = m = and n = I = 1, the above bracket returns Darcy's law for p = Aq in 
the more general form 

dp{q, t) d ( , . f , , , . d 6E 



dt dq 



(^p{q,t) I K,^,{q,q') p,{q',t) dq'^ (2.14) 



which recovers (j2.8p when Ki^i{q,q') = 6{q — q'). It is interesting to notice that, in the multi- 
dimensional treatment, the existence of the kernel i^i,i allows for a mobility matrix, rather than 
just a single-valued mobility functional. Indeed, the quantity is always allowed to involve 
a symmetric positive matrix Ki^i, as it can be seen by rewriting the Taylor expansion ()2.13p in 
the form Er(q, p, q', p') = ^p®" (g) p'®'"_| i^^ ,„(q, q'), where the symbol _l denotes tensor 
contraction and Kn,m is a covariant tensor field of order n + m. Thus Darcy's law becomes, for 
example, 

I = d,v V) ^) =div(p(A,V)f) (2.15) 

where we have introduced the matrix Ki^i so that -ftri^i(q, q') = /t'l^i (5(q — q') and we have 
defined the mobility matrix pp := ppKi^i. Physically, this corresponds to the case of isotropic 
particles moving through a non-isotropic medium. 

The remainder of the paper concentrates on the case when the sharp operator ( • )" is defined 
in terms of a delta-function kernel and the mobility pp is a real-valued functional of p. 



3 Dissipative dynamics with anisotropic interaction 

3.1 Purely dissipative Vlasov equation using GHK bracket 

Following GHK, we introduce a particle distribution which depends not only on the position 
and momentum coordinates q and p, but also on an extra coordinate g associated with, e.g., 
orientation or some other order parameter. The coordinate g belongs to the dual q* of the 
Lie algebra g of the order parameter group G, which for anisotropic interactions is ^ = S0{3). 
For an introduction to this approach in the theory of complex fluids, see |Ho2002) . In what 
follows, we shall formulate the problem in the general context and analyze the case of rotations 
separately. In the non-dissipative case, the Vlasov equation is written in terms of a Poisson 
bracket, which is the direct sum of the canonical ((7, j>)-bracket and the Lie- Poisson bracket on 
the dual Lie algebra g*. Explicitly, this Poisson bracket is written as (|1.3p . Like in the isotropic 
case, the non-dissipative Vlasov equation now becomes 

dt 

Analogously to the isotropic case, the above Vlasov equation is determined by the Lie-algebra 
action h-f := {f,h}^ of phase-space functions h on phase-space distributions /. The correspond- 
ing dual operation {*) is given by {f -k k,h) = {k, — h ■ f) = ( k, {h, f}^ ) = ( {/, k}^ , h). Here 
we have denoted by (•, •) the pairing between functions and distributions on phase space and 
the equality {k, {h,f}^ ) = {{f,k}-^ , h) follows from Leibniz rule for Poisson brackets. Thus, 
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= {f,k}^. Upon applying the same arguments as in the previous Section [2] and making use 
of the general formulation of the dissipative bracket (jl.Sp . we find the purely dissipative Vlasov 
equation in double-bracket form (anisotropic DBV equation): 



where we have dropped the sharp symbol. This equation has exactly the same form as in ()1.2p . 
but now one replaces the canonical Poisson bracket {• , •} with the direct sum Poisson bracket 
{■ 1 'li in (|l-3p . This formulation can now be used to derive the double-bracket dissipative version 
of the Vlasov equation for particles undergoing anisotropic interaction. 

3.2 Dissipative moment dynamics in ID 

An extension of the ID moment dynamics hierarchy to include orientation dependence, may 
be obtained by following the same steps as in the previous section, beginning by introducing 
the moment quantities An{q,g) in (II. 4p . One may find the entire hierarchy of equations for 
these moment quantities and then integrate over g in order to find the equations for the mass 
density p{q) := j Ao{q, g) dg and the continuum charge density G{q) = J g Ao{q, g) dg. Without 
the integration over g, such an approach would yield the Smoluchowski approximation for the 
density AQ{q,g), usually denoted by p{q,g)- This approach is developed in Sec. HJ where the 
dynamics of p{q, g) is presented explicitly. 

In this section, we extend the Kupershmidt-Manin approach as in GHK to generate the 
dynamics of moments with respect to both momentum p and charge g. The main complication 
is that the Lie algebras of physical interest (such as so (3)) are not one-dimensional and in general 
are not commutative. Thus, in the general case one often introduces a multi-index notation as 
in |Kul987l IGiHoTr2005j . Rather, we introduce tensor powers g" := g®'^ = g ® ■ ■ ■ ® g. Then, 
the moments are expressed as 



Note that all the integrals with respect to dp and dg are taken in finite-dimensional spaces and are 
presumed to exist according to Assumptions 1 2 . 1] and 12 . 2 1 This multi-dimensional treatment leads 
to cumbersome calculations |GiHoTr2008j . For the purposes of this section, we are primarily 
interested in the equations for p and G; so we shall restrict our considerations to moments 
of the form An^ = J p"^ f{q,p, g) dpdg and An^i = J p"^ g f{q,p, g) dpdg. We write the single 
particle Hamiltonian as h = 6H/6f = p"'g'^ -1 5H/6An^u ='■ P^g'^ I3n,u{q)- Here _l denotes 
tensor contraction between lower and upper indexes. We restrict our considerations to the case 
u = 0,1. Consequently, we employ the following assumption. 

Assumption 3.1 The single-particle Hamiltonian h = 5H/5f is linear in g and can be ex- 
pressed in ID as h{q,p,g) = p^ ipn{q) +p"'{g, V'nl?))? where ipniq) G M is a real scalar function 
and ip,^{q) € M®g is a real Lie-algebra-valued function. This assumption will be used throughout 
the rest of the paper. 

3.2.1 Dual Lie algebra action 

The action of /3n,i/ on / is defined as (3n,u • / = {p"g'^ -i f^n,i^, f}i (no sum), where _l denotes 
tensor contraction between lower and upper indexes. The dual of this action is denoted by 






(3.2) 
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{*n,u)- It may be computed analogously to (|2.5|) and found to be 



jj P'^d'y, k\^P^9 
/ 5^+'" J ad^,^^ dg + J g" (g, 









5g 



dg 
dg. 



dg 



(3.3) 



Here, k = g'^ ^ ara,a{<i) (with a = 0, 1) and we have used the definition of the moments 
An{q-, g) in (ll.4p . Moreover, the notation ad*^ ^An^a+u above involves a tensor index contraction 
(arising from the contraction symbol in the second line of (j3.3p ) between the upper and lower 
Lie algebra indexes. For example, in the simplest case one has ad*gj^o,i = (^o,i)a 9g(ao,i)"- 

Once the dual of the Lie algebra action is completely characterized, we can now write the 
double bracket structure for the moments An^u- This is done in the following discussion. 

3.2.2 Evolution equations 

Having characterized the dual Lie algebra action in (|3.3p . we may write the evolution equation 
for an arbitrary functional F in terms of the dissipative bracket as follows: 



F = {{F,E}} 




(3.4) 



where the pairing ((•,•)) is given by contraction of all upper and lower tensor indexes and inte- 
gration over the spatial coordinate q. In the one dimensional case under consideration, we have 
explicitly / •kn,i' SF/6f = (/ -kn,^ 5F /5f^ . . e*^ (g) . . . (g) e*" , with |e*| being the basis in the 
dual of Lie algebra. 

We now consider the lowest order by setting n = 0, 1 and excluding the case (?i, v) = (1, 1) (we 
have 1/ = 0, 1 as before). This is exactly the same truncation that is involved, on the Hamiltonian 
side, in the Vlasov moment equations of chromohydrodynamics |GiHoKul982l IGiHoKul983j . 
Thus, upon denoting k = 6F/6f, we evaluate 

/ *0,0 k = ad^^^^ Am-l, a = ^^Oimfi + ^^am,l "^^-1, 1 





dA^ 


[9(9, 


dg 



d9 = 3d*a^o Am-i, 1 + ad*a^ ^ Am,i 



f *i,o k — ad% Am, a 



ad 



a, 



(3.5) 



where k = p™ g'^ _l am,a (<?) and we have truncated to consider only moments of the type Amfi 
and A„i,i- Note that the notation for the Kupershmidt-Manin ad* operator differs from the 
notation for the coadjoint operator ad* associated to the Lie algebra g. At this point, we 
consider the simplest case in which m = 0, so that the first two relations of ()3.5p are /*o,o k = 
and / *o,i k = ad*^ ^ Aq^i, while the third of ()3.5p becomes 



/★i,o k = ad^Qo Ao^o + ad^g^^ Aq^i = Aq^q dgOofi + (^o,i)a9g(ao,i)" 
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The equation for the evolution of -F = Aq x is found from the bracket l\'6A\) in the truncated form 



/^fo, \- { fJ-o,! , — ) , Mo, ^ {Ao,i, 



oq \ c)q / J \ oq \ oq 

^ad^ci/^o,!)", ad;^^^o,i 



where the prime in the summation symbol stands for n ^ v and we have expanded 5E/5f = 
/3o,o + (5,/3o,i> = / 5AQfl + {g,5E / 5Aq^i) . Also, we use the notation = / p"^ fi[f] dp dg . 
The sharp operator in the above formulas involves an appropriate metric: all the discussions 
from section [23] involved in the derivation of Darcy's law also hold in this context. 

We now simplify the notation by defining the following dynamic quantities: p = Aq^q and 
G = Aq^i. Likewise, we define the mobilities as: pp = /io,o and pc = fJ-o.i- In terms of these 
quantities, we obtain the following result. 

Theorem 3.2 The ID moment equations for p and G are given by 
dp d ( ( d 5E / d 5E\ 

Tt = dqy[!''d-q-p^Y''^d-qJG^^ ' ^^'^^ 




Equations in this family (called Geometric Order Parameter equations) were derived via a dif- 
ferent approach in |HoPu2005| IHoPu2006t IHoPu2007] . These equations were also derived by 
using a multi- index approach in |HoPuTr2008j . 

Remark 3.3 (Higher dimensions) Gibbons' Theorem \2.3\ has recently been extended to ac- 
count for more geometric structure in I GiHo Tr2008^ . where the authors also provide a general 
setting for the multi- dimensional treatment. The higher dimensional treatment of the moment 
Poisson structure provides a systematic machinery that enables us to write the double bracket 
moment equations in the full three dimensional case. In 3D, the equations ()3.6p and (j3.7p become 

dp ( ( I bE 

-=div yp (/^pV- + (^^c,V- 

f ^div(c(.,vf.(..,V§;jJ.ad*^^^^^^^,C, (3.8) 
For the seek of simplicity, this paper considers only the one dimensional case. 




3.3 Singular solutions 

Equations ()3.6p and ()3.7p admit singular solutions. This means that the trajectory of a single 
fluid particle is a solution of the problem and all the microscopic information about the particles 
is preserved. We shall prove the following. 
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Theorem 3.4 Equations 6\) and ^3. 7p admit solutions of the form 



p{q,t) = Wp(t) 5{q - Qp{t)) 



G{q,t)=WG{t) 5{q-QG{t)) 



(3.9) 



where liip = 0, wq 



ad*jj wg and 



d 6E 



fJ-G, 



d_5E_ 



G 



d 6E 



f^G 



d 6E 



q=Qg 



A proof of this theorem is easily obtained by taking the pairing of equations (j3.6p and (j3.7p 
with a couple of test functions {4^ p{q) ■, (I)g{<1)) ■, integrating by parts and then substituting the 
singular solution anstatz. Again, we shall assume that all the necessary integrals exist since 
test functions {(j)p{q), 4>g{q)) ai'e smooth, supposed to decay to zero sufficiently fast for large q 
if necessary and {p, G) together with the first derivatives are in the generalized function space. 

A similar result applies for the Geometric Order Parameter (GOP) equations investigated 
in [HoPu2005[ rHoPu2006[ lHoPu2007| . 



3.4 An illustrative application: a filament of rod-like particles 



In this case G = m(x), adyW = v x w and ad*w = — v x w, and the Lie algebra pairing is 
represented by the dot product of vectors in M'^. Therefore the equations are 



dp d 



d 6E 



d 6E^ 



dt dx\^ \ dx 5p ~'~ dx Sm.^ 



dm 



d_ 

dx 



m 



d 6E 
^^d-xJ-p^^^ 



d SE\ 
dx 6m J 



+ m X 1/1 



SE 
6m 



(3.10) 
(3.11) 



Note that equations for density p and orientation m have conservative parts (coming from the 
divergence of a flux). In addition, when fi^ = am for a constant a, the orientation m has 
precisely the dissipation term m x m x 6E/6m introduced by Gilbert [Gilbert 1955j . Thus, we 
have derived the Gilbert dissipation term at the macroscopic level, starting from double-bracket 
dissipative terms in the kinetic theory description. As far as we know, this is the first time that 
the Gilbert dissipation term has been derived from a kinetic theory model. 



4 Smoluchowski approach 

We shall now turn our attention to the Smoluchowski approach for the description of the in- 
teraction of anisotropic particles. Usually, these particles are assumed to be rod-like, so their 
orientation can be described by a point on a two-dimensional sphere [DoEdl988| . In this 
picture, one defines a distribution p{q, fi, t) on the extended configuration space x S^, whose 
dynamics is usually regulated by the Smoluchowski equation 



^ = divq (^pht{p) + divn (^p ^lo{p) (^-l) 



where p,T and po are respectively the mobilities in the translational and orientational part and 
the rotational gradient operator Vn is expressed as Vn = fi x d/dfl. Upon combining the 
previous derivation of Darcy's law with the definition of the two different mobilities, one easily 
recognizes that the above equations can be obtained by following the same steps as in Darcy's 
law. Two different mobilities pT and po involving different smoothing kernels are allowed by 
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the particular Cartesian product structure of the configuration space M'^ x S'^. The reference 
Vlasov equation for the above Smoluchowski equation is written as 



dt \ { 6f } -p.jja I of J 



T*(M3x52) 



This principle could also be used to define two different mobihties /x^ and /x^ in equation (|3.1ip , 
so that /x^ would only appear in the expression of Darcy's velocity, while /x^ would appear in 
the Gilbert dissipation term involving the double cross product. For the sake of simplicity, we 
shall not deal with such a distinction between these two types of mobilities. 

In what follows, we shall consider particles of arbitrary shape, for which one needs the full 
group 5*0 (3) to define their orientation. We shall work with the corresponding Lie algebra 
so (3) to conform to our theory. This construction particularly suits the dynamics of ferro- 
magnetic particles. Following the Smoluchowski approach, moments are defined as An{q,g) ■= 
f P"' /{q^Pt 9) dp- As in the Kupershmidt-Manin approach, these moments are dual to /3n{q,g), 
which are introduced by expanding the Hamiltonian function h{q,p, g) as h{q,p, g) = /3n{q, g)- 
Substitution of these forms yields a Lie algebra bracket for the quantities /3„ given by |/3„, am\i - 
I/^njOm] + {g, [/^njC^m])' where prime denotes partial derivative with respect to g and we have 
used the same notation as in Section 12.21 A substantial difficulty arises in this moment hier- 
archy because the above Lie algebra structure does not allow for subalgebras, as happens for 
the moments treated in the previous sections. Thus, there is no rigorous way of truncating 
the hierarchy, which instead needs to be truncated by ad hoc physical arguments. However, in 
what follows we notice that Smoluchowski equation (for 50(3) orientations, or ferromagnetic 
particles) can be obtained by taking the zero-th order moment AQ{q,g) of the oriented DBV 
equation. 

4.1 A cold-plasma closure and its dissipative dynamics 

We shall apply the double bracket approach to the following equations for the fluid momentum 
i^{q) = / P f{q-,P-,g)^P^g and the probability distribution p{q,g) = f f{q,P,g)dp = Ao{q,g) in 
the orientation 5 € 0*: 



dt 6m dq dq 5m J ^ dq 5p ^ ^ dt dq \ 6m J ^' 

These equations arise from the Vlasov cold-plasma solution 

m{q) 



dp d_m 

dg ' dg 6p 



f{q,p,9) = p{Q,g) Hp 



fp{q,g) dg 



and their Hamiltonian structure is directly inherited from the GHK- Vlasov bracket jGiHoKul983j 
{F, H} = J f ^6F /6f , 6H /6f^-^dqdpdg ,hy direct substitution of the chain rule formula = 
6F/6p + p 6F/6m for generic functionals F and H. 

In this picture, the Lie algebra action of the variable {^fiq, g),u{q)) S C°°{Q x g*) x X{Q) 
is given by (93, u) ■ f = {/, ip + pu}^ and its dual action is given by 

(^f-kk, {p,u)^ := -{k, {lp,u) ■ f^j = (^(^j {f,k}idp, j p{f,k}idpdg^ , {p,u)^ , 

where the last pairing (• , •) : {C°^{Qxq*)xX{Q))x{C°°{Qxq*)xX{Q))* ^ M between fluid and 
particle variables is defined as the sum ((p, m), {(p, n)) = / p{q, g) <p{q, g) dqdg + f m{q) u{q) dq. 
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More explicitly, the star operator is defined for k{q,p,g) = il'{(l-,g) + pv{q) as 



( d{pv) 



dp dip 



dm dv f dip 
, V— + 2m— + p — dg 
dq dq J dq 



(4.3) 



At this point, we are ready to introduce the dissipative bracket by 



F = {{F,E}} 




dpp d 6E 
dg ' dg dp 



'^^^ dq dp ^ 



dp d 6F' 
dg ' dg dp 



^ dq dp 



(4.4) 



where pp = j pdp. Here we have restricted to consider only functionals of the probabihty 
distribution p{q,g), so that 6E/6f = 6E/6p and analogously for F (no inertial effects). By 
using this evolution equation for an arbitrary functional F, the rate of change for zero-th moment 
p{q, g) = I f{Q, P, a) dp is found to be 



di 

dt 



d_ 
dq 



d 5E ' 
^''d-q-p'^'. 



5E\ 



(4.5) 



where we have introduced the notation {•, - jg* for the Lie-Poisson bracket {(/9, V'jg. := (5, [dg ^ , dg 
The novel feature of equation ()4.5p is the integral appearing in the expression of Darcy's velocity 
(first term of the right hand side). While ordinary Smoluchowski equations treat the transla- 
tional and orientational coordinates on the same level, equation (j4.5p treats these coordinates 
from two very different perspectives, while retaining the Yang-Mills charge nature of the orien- 
tation coordinate g. This reflects in the appearance of the integral over the orientation g that 
produces a velocity vector field defined on the manifold Q. Such an integral is not present in 
conventional Smoluchowski approaches. 

Remark 4.1 (Higher dimensions) Analogously to remark \3. 3|, we can also write the three 
dimensional version of equation (j4.5p . This reads as 



dp 

m 

4.2 Singular solutions 



div \^p 



A^P V— d^ 




Both conventional Smoluchowski equation (j4.ip and equation ()4.5p possess delta-function solu- 
tions of the type p{q,g) = T S{q — Q{t)) 6{g — G{t)), where T is a constant. This is the usual 
Klimontovich single particle solution. It is an open question whether such solutions emerge 
spontaneously. Another class of singular solutions is the following 

p{q,g)=w{g,t)6{q-Q{t)). (4.6) 

This is a point particle solution that supports a kinetic equation for its own probability distri- 
bution w{g, t) in the space Den(g*) of densities in the orientational degree of freedom 5 S 0*. In 
this case, the dynamical equations are 



Q 



d_5£ 
dq 6p 



q=Q 



dw 
'dt 




This is a coupled system for the ODE regulating the particle trajectory Q{t) and the PDE 
governing the orientational probability w{g, t). To our knowledge, this class of singular solutions 
is new. 
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5 Summary and outlook 

The double-bracket Vlasov moment dynamics discussed here has provided an alternative to both 
the variational-geometric approach of [HoPu2007j and the Smoluchowski approach reviewed in 
|Co2005j . These are early days in this study of the benefits afforded by the double-bracket 
approach to Vlasov moment dynamics. However, the derivations of the Darcy law in (12. Sh and 
the Gilbert dissipation term in (|3.1ip by this approach lends hope that this direction will provide 
the systematic derivations needed for modern technology of macroscopic models for microscopic 
processes involving interactions of particles that depend on their relative orientations. Although 
some of these formulas may look daunting, they possess an internal consistency and systematic 
derivation that we believe is worth pursuing further. Our next steps will be the following: 1) 
perform the analysis of the mobility functionals in kinetic space //[/] as well as the mobilities for 
each particular geometric quantity /Xp, ^g, etc.; 2) determine the conditions for the emergence 
of weak solutions (singularities) in the macroscopic (averaged) equations. 
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